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Breaking of N = 8 magicity in J Be 
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Structure of 13 Be was investigated with antisymmetrized molecular dynamics. The variation after 
spin and parity projections was performed. An unnatural parity 1/2" state was suggested to be 
lower than 5/2 + state indicating that vanishing of the N = 8 magic number occurs in 1 Be. A 
low-lying 3/2 + state with a 2hui configuration was also suggested. Developed cluster structures 
were found in the intruder states. Lowering mechanism of the intruder states was discussed in terms 
of molecular orbitals around a 2a core. 
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I. INTRODUCTION 

One of the exotic phenomena discovered in unsta- 
ble nuclei is vanishing of neutron magic numbers. In 
neutron-rich Be isotopes, the breaking of N = 8 magic- 
ity in 11 Be has been known from the abnormal spin and 
parity 1/2+ of the ground state. The vanishing of the 
N = 8 magic number in 12 Be has been suggested by slow 
j3 decay to 12 B vA and it has been supported by various 
experiments 040] ■ The N = 8 shell breaking has been 
suggested also in n Li by experimental and theoretical 
works 01 ■ Even though 12 Be and n Li are neighboring nu- 
clei, they have different characters and the shell breaking 
mechanism is not the same between two nuclei. One of 
the remarkable features of 12 Be different from n Li is de- 
veloped cluster structure, which plays an important role 
in the shell breaking in Be isotopes. 

Cluster structures of Be isotopes have been intensively 
investigated in many theoretical works [8H29J] . Low- lying 
states in neutron-rich Be were described successfully with 
cluster models and molecular-orbital models assuming 
two a clusters and surrounding neutrons. The forma- 
tion of a 2a core in Be isotopes was confirmed by the 
author and her collaborators with a method of antisym- 
metrized molecular dynamics which does not rely on a 
priori assumption of any clusters J14l4l9{ . These works 
revealed that 2a structures are favored in neutron-rich 
Be isotopes where valence neutrons around the 2a core 
play important roles. To understand cluster features of 
low-lying states of neutron-rich Be, a molecular orbital 
picture is helpful |8T-[l3|. In the picture, molecular or- 
bitals arc formed by a linear combination of p orbits 
around two a clusters, and valence neutrons occupy the 
molecular orbitals. A longitudinal positive-parity orbital 
is called "er^-orbital" . Since the tTi/ 2 -orbital has two 
nodes along the a-a direction, it gains kinetic energy as 
the 2a cluster develops. Note that the cr^-orbital cor- 
responds to the "1/2 [220]" orbit in the Nilsson model 
(a deformed shell model) [301 an d it originates in the sd- 
orbit in the spherical shell model limit. In the molecular 
orbital picture, the ground states of 11 Be and 12 Be are 
regarded as the configurations with one and two neu- 
trons in the c^-orbital, respectively, which correspond 
to the intruder configurations in terms of the spherical 
shell model. In other words, lowering mechanism of the 
intruder states in Be isotopes can be understood by the 



the energy gain of the er^-orbital in the developed 2a 
cluster structure. 

It is then a challenging issue to investigate cluster fea- 
tures and shell evolution in further neutron-rich Be iso- 
topes near the drip line. Indeed, many experiments have 
been performed to observe energy spectra of 13 Be (see 
Refs. [3l|,|32[ and references therein). Since 13 Be is an un- 
bound nucleus, all energy levels are resonance states. A 
resonance about 2 MeV above the neutron-decay thresh- 
old has been observed in several experiments and it has 
been assigned to a 5/2+ state J3ll43a |. However, other 
levels of 13 Be have not been confirmed yet. The position 
of a 1/2" state relative to the 5/2 + state is, in particular, 
a key problem concerning the breaking of N = 8 magicity. 
Recently, energy spectra have been measured with invari- 
ant mass measurements of 12 Be+n by Simon et al. |3l| 
and Kondo et al. [32J. In the former work, they argued 
that an observed low-energy peak just above the 12 Be+n 
threshold is described by s-wave virtual state contribu- 
tion. On the other hand, the latter work reported the 
existence of a p-wave resonance in the low-energy peak 
below the 5/2 + state and tentatively assigned it to a 
1/2" state suggesting the breaking of N = 8 mag icity 
in 13 Be. However, a controversial claim was given 



;icity 



state is still under 



and the energy position of the 1/2 
discussion. 

If the inversion of single-particle levels does not occur 
in 13 Be, the unnatural-parity 1/2" state should be higher 
than natural-parity l/2 + and 5/2 + states. The question 
is whether or not the inversion occurs and the l/2~ state 
comes down to the low-energy region in 13 Be. Provided 
that the inversion occurs, cluster structure may play an 
important role in the intruder 1/2" state as well as n Be 
and 12 Be. Some theoretical calculations sugges ted possi- 
bility of a low-lying 1/2" state in 13 Be[32|, |38|. For in- 
stance, a shell model calculation for 13 Bc using the SFO 
interaction [40fl , which is adjusted to reproduce the par- 
ity inversion of 11 Be, gives the 1/2" ground state in 13 Be 
[32j |. In microscopic cluster model calculations, there is 
only few applications to 13 Be [39[ because of the large 
number of valence neutrons. 

Our aim is to investigate structure of 13 Be. A method 
of energy variation after spin-parity projection (VAP) in 
the framework of antisymmetrized molecular dynamics 
(AMD) 0,111, El| is applied to 13 Bc. The AMD+VAP 
method has been already applied for studying 10 Be[L7|- 



n Be[l8|, and 12 Be[19[ and described successfully proper- 
ties of ground and excited states. In the energy levels of 
13 Be calculated by using the interaction that reproduces 
the parity inversion of 11 Be, we will suggest a low-lying 
l/2~ state below the 5/2 + state. We show the coexis- 
tence of Ofuu, Yfau), and 2Tiuj states in a low-energy region 
of 13 Be and suggest the breaking of N — 8 magicity. We 
discuss cluster structures of 13 Be and indicate that the 
calculated states of 13 Be can be classified by molecular 
orbital configurations. The neutron shell breaking is dis- 
cussed systematically in a chain of Be isotopes focusing 
on cluster aspect. 

This paper is organized as follows. In the next section, 
the formulation of the present calculation is explained. 
The results are shown in section IIIII and discussions are 
given in section ITVl Finally, in section [VI a summary and 
an outlook are given. 



II. FORMULATION 

We describe 13 Be with AMD wave functions by apply- 
ing the VAP method. For the AMD+VAP method, the 
readers are referred to Rcfs. [17H191 |42| . The method 
is basically the same as those applied to 10 Be, n Be, 
and 12 Be. A difference in calculation procedures from 
Rcfs. |17H19| is that we do not adopt an artificial bar- 
rier potential in the present calculation which was used 
in the previous works. To see effects of quasi-bound fea- 
tures of the last valence neutron on energy spectra, we 
also adopt 12 Be+n wave functions by using the 12 Be core 
wave functions obtained with AMD+VAP for 12 Be. 



A. AMD wave functions 

An AMD wave function is given by a Slater determi- 
nant of Gaussian wave packets; 
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where the «th single-particle wave function is written by 
a product of spatial(^), intrinsic spin(x) and isospin(r) 
wave functions as, 
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<hti an d Xi are spatial and spin functions, and Tj is 
the isospin function fixed to be up (proton) or down 
(neutron). Accordingly, an AMD wave function is 
expressed by a set of variational parameters, Z = 
{Xi, X 2 , • • • , Xa,£i,£2, • • • , £a}- The width parameter 
v is chosen to be v = 0.17 fm -2 which is the same value 
as that used for 12 Bc in Ref. (Tal. 



B. Variation after projection method 

Energy variation after spin and parity projections 
(VAP) in the AMD model space is per formed as is done in 
the previous studies of Be isotopes |17H13 | . For the lowest 
J* state, the parameters X, and £i(i = 1 ~ A) are varied 
to minimize the energy expectation value of the Hamil- 
tonian, ($|iJ|<t>)/(<f>|$), with respect to the spin-parity 
eigen wave function projected from an AMD wave func- 
tion; $ = P//^-<I>amd(Z). Here, P(f K is the spin-parity 
projection operator. The energy variation is performed 
with a frictional cooling method [16[ . Then the optimum 
AMD wave function <E>AMD(Z j7r ), which approximately 
describes the intrinsic wave function for a J 7 * state, is 
obtained. For each </ 7r = J^ a , the optimum parame- 
ters Z jQ7r ° are obtained. After the VAP procedure, final 
wave functions are calculated by superposing the spin- 
parity eigen wave functions projected from all the AMD 
wave functions <I>AMD(Z jQ7rQ ) obtained by VAP for vari- 
ous J^ a states. Namely, the final wave functions for the 
J™ states are expressed as, 
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where the coefficients c J (K,J a ,7T a ) are determined by 
diagonalization of norm and Hamiltonian matrices. 



C. 12 Be+n model 

As shown later, the results obtained with AMD+VAP 
show that low-lying states of 13 Be can be interpreted as 
12 Be+n, where the 12 Be core is the intrinsic state of the 
12 Be(0^") having an intruder 2froj configuration or that 
of 12 Be(0^) with a normal Ohm configuration. In reality, 
13 Be is an unbound nucleus and all states are resonances 
above the 12 Be+n threshold. In such a case, asymp- 
totic behavior of the valence neutron wave function in the 
outer region can be important, in particular, for energy 
position of low angular-momentum states. However, the 
AMD method is not suitable to describe detailed behav- 
iors of asymptotic regions, because a system is expressed 
by a Slater determinant of Gaussians and is treated in a 
bound state approximation in the AMD model. To see 
effects from spatial extension of the last neutron on en- 
ergies of resonances, we also calculate energy levels in a 
12 Be+n model described below. 

In the 12 Be+n model, wc first apply AMD+VAP to 
12 Be to obtain intrinsic wave functions $( 12 Be : 2huj) 
of 12 Bc(0f ) and $( 12 Bc : Oftw) of 12 Bc(0^) 0. Each 
intrinsic wave function $( 12 Bc : a) (a = Ohui , 2hui) is 
expressed by a AMD wave function and it is written by a 
Slater determinant. Then we add one neutron to the core 
wave functions $( 12 Be : 2Huj) and $( 12 Be : 0hu>). The 
additional neutron wave function in the 12 Be+n system is 
described by a Gaussian wave packet located at a position 
X relative to the core. A 13 Bc wave function for a J v 



state is described by a linear combination of 12 Be+n wave 
functions with various positions X as 
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where VV|-(+) ( x ) i s a spin-up (spin-down) neutron wave 
function with a Gaussian form with the width v param- 
eter localized at X. 
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The 12 Be wave function is shifted by — X^/13 to take 
into account the recoil effect from the last neutron. For 
the total wave function, the antisymmetrization and the 
spin-parity projection are performed as well as the su- 
perposition of basis wave functions. 

In the present calculation, the intrinsic states $( 12 Be : 
a) are deformed and their orientations are chosen to sat- 
isfy (.t 2 ) < (y 2 ) < (z 2 ) and (xy) = (yz) = (zx) = 0. 
For the position X& of the last neutron Gaussian wave 
function, grid points in the |x| < 5 fm and \z\ < 5 fm 
region on the y — plane are taken. The grid spacing is 
chosen to be 1 fm. 

When an intrinsic wave function of the 12 Be core is 
axial symmetric, X& on the y plane is enough to take 
into account coupling of the last neutron with all rota- 
tional band members of 12 Be constructed from the in- 
trinsic wave function. Strictly speaking, the present 12 Be 
wave functions are not axial symmetric, however, X& is 
restricted only on the y = plane to save numerical cost 
in the present calculation. 

With the 12 Be+n model, we calculate 13 Be energy 
spectra. Comparing the 12 Be+n model calculation with 
the AMD+VAP one, we will discuss, in particular, how 
the energy spectra can be modified from the AMD+VAP 
results by improving the last neutron wave function. 



III. RESULTS 

A. Effective interaction 

We used the same effective nuclear interaction as that 
used in Rcfs. [H Ql]. It is the MV1 force JH for 
the central force supplemented by a two-body spin-orbit 
force with the two-range Gaussian form same as that 
in the G3RS force |4J|. The Coulomb force is approx- 
imated using a seven-range Gaussian form. We adopted 
the interaction parameters that are used for 11 Be and 
12 Bc [H, [H- Namely, the Majorana, Bartlett, and 
Heisenberg parameters in the MV1 force are m = 0.65, 
b = 0, and h = 0, respectively, and the spin-orbit 



strengths are taken to be m = —ujj = 3700 MeV. We 
denote this parameterization by the set (1) interaction. 
To see the interaction dependence of the theoretical re- 
sults, we also used the other parametrization (2) with 
weaker spin-orbit forces, uj = —un = 2500 MeV. 

The energy levels of the excited states of 10 Be, 11 Be, 
and 12 Be are reproduced well by the AMD+VAP cal- 
culations with the set (1) interaction. Particularly, the 
breaking of the N = 8 magicity in 11 Be and 12 Be are suc- 
cessfully described with the set (1) interaction (l8l KL9|. 
In this paper, we mainly discuss the results calculated 
with the set (1). 



B. Energy levels 
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FIG. 1: Energy levels of (upper) 13 Be and (lower) 11 Be. Exci- 
tation energies are shown. The calculated levels are obtained 
with AMD+VAP by using the set (1) and (2) interactions. 
Results for "Be obtained with AMD+VAP are taken from 
[la |. The experimental data of 13 Be are those in Ref. |33 |. 



> 






1/2 + 




3/2" j 




5/2" :' 


- 


5/2* 

3/2+ • 


A 


" 


ly 3/2" — 

H5/2" 

t 3 ^- 5/2 + 

3/2 + l/2+ 5/2 + _ 


. 


1 5/2 + 

1/2" 




1/2" 


VAP(l) 


12 Be+n exp(a) exp(b) 



FIG. 2: Energy levels of 13 Be calculated with AMD+VAP 
and those with the 1 Be+n model by using the set (1) inter- 
action. Energies relative to the 12 Be-n threshold are shown. 
The calculated energies are measured from the theoretical 
12 Be-7i threshold energy, —61.3 MeV. Experimental energies 
are those taken from (a) Ref. [32|] and (b) Ref. [31| ] . 



We applied AMD+VAP to 13 Be and calculated J n 
states up to J = 5/2. The calculated energy levels of 
13 Be are shown in Fig. Q] compared with experimental 
energy levels reported in Ref. [32[. Ener gy l evels of n Bc 
are also shown. As mentioned in Ref. [lq . the set (1) 
interaction reproduces the unnatural parity ground state 
1/2+ in 11 Be, while the set (2) interaction fails to de- 
scribe the parity inversion. 

In the result with the set (1) interaction, energy levels 
of 13 Be are found to be out of the normal ordering. 1/2" 
and 3/2 + states almost degenerate at the lowest energy 
and a 5/2+ state exists above them. As discussed later, 
the 1/2" and 3/2 + states have dominantly Vfrw and 2huj 
excited configurations, respectively, while the 5/2 + state 
is described by a normal Ohw configuration. The appear- 
ance of these intruder states in such the low-energy region 
suggests the breaking of N = 8 magicity in 13 Be as well 
as n Be. 

In contrast, in the result with the set (2) interaction, 
the 5/2+ state is the lowest and the l/2~ and 3/2+ states 
are higher than it as naively expected from the spherical 
shell model. 

Although all states are resonances above the 12 Be+n 
threshold in a real 13 Be system, they are treated in a 
bound state approximation in the present AMD+VAP 
calculation. In principle, a resonance energy might de- 
crease when the last neutron occupies a low angular- 
momentum orbit because of an extended neutron wave 
function in an outer region. It means that asymptotic 
behavior of the valence neutron should be taken into ac- 
count carefully for more detailed discussion of energy lev- 
els. To sec how the level ordering is affected by improving 
wave functions for the last neutron, we calculated energy 



levels of 13 Bc also in the 12 Be+n model. 

The 13 Bc energy spectra calculated with the 12 Be+n 
model by using the set (1) interaction are shown in Fig. [3] 
The calculated energies are measured from the theoret- 
ical value —61.3 MeV of the 12 Be-n threshold energy, 
which is evaluated by diagonalizing spin-parity eigen 
states projected from $( 12 Be : 2%w) and $( 12 Be : OTiw). 
For all states, the 12 Be+n model calculation gives about 
a few MeV lower energies than those obtained with 
AMD+VAP because of improving neutron wave func- 
tions. The level ordering somehow changes from the VAP 
results, for instance, the 3/2+ state shifts up while the 
l/2 + state comes down relatively. However, the 12 Be+n 
calculation shows again the feature of the neutron magic 
number breaking that various spin and parity states de- 
generate in the low energy region. In particular, it should 
be pointed that the 1/2" state is the lowest consistently 
with the experimental report by Kondo el al. \32\ . 

In experimental and theoretical studies of 13 Bc, the 
l/2 + state has been suggested to be a virtual state and 
contribute to the spectra near the 12 Be+n threshold en- 
ergy [3l|, [33, [39| ■ In the present calculation, model space 
is not enough to describe a virtual state. That may be the 
reason why the l/2 + state is still higher than the 1/2" 
and 5/2 + states even in the 12 Be+n model calculation. 



IV. DISCUSSIONS 

In the present AMD+VAP calculation, any clusters 
are not assumed in the model. Nevertheless, the results 
suggest that cluster structures appear in 13 Be as well as 
other Be isotopes. Moreover, the 13 Be states obtained 
with AMD+VAP can be associated with 12 Be+rt states, 
and their structures correspond to those of the 12 Be+n 
model calculation. Since each intrinsic wave function ob- 
tained by AMD+VAP is expressed with a single Slater 
determinant, AMD+VAP wave functions are useful to 
analyze intrinsic structures. We here discuss cluster fea- 
tures of 13 Be by investigating AMD+VAP wave func- 
tions while focusing on the 2a cluster and valence neu- 
tron configurations. We also discuss systematics of low- 
lying states in neutron-rich Be isotopes in the molecular 
orbital picture. 



A. Intrinsic structures and clustering 

As seen in density distributions of intrinsic wave 
functions shown in Fig. [3] a 2a core is formed in 
13 Be. In particular, developed 2a cluster structures with 
large deformations are found in negative-parity states, 
13 Be(l/2"), 13 Be(3/2"), and 13 Be(5/2"). In positive- 
parity states, 13 Be(l/2 + ) and 13 Be(5/2 + ) show relatively 
weaker cluster structures with smaller deformations while 
13 Be(3/2 + ) shows a remarkable 2a structure similarly to 
the negative-parity states. 
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FIG. 3: (Color on-line) Density distributions of the intrinsic 
states for l/2~, 3/2~, 5/2", 1/2+ , 3/2+ , and 5/2+ states 
in 13 Be calculated with AMD+VAP using the set (1) inter- 
action. The orientation of an intrinsic state is chosen so as 
to satisfy (a; 2 ) < (y 2 ) < (z 2 ) and (xy) = (yz) = (zx) = 0. 
The horizontal and vertical axes are set to the z and y axes, 
respectively. Densities are integrated with respect to the x 
axis. Distributions of matter, proton and neutron densities 
are shown left, middle, and right, respectively. 



Since an intrinsic wave function for each state is given 
by a Slater determinant in the AMD+VAP calculation, 
we can analyze single-particle wave functions and discuss 
neutron configurations as well as single-particle energies. 
Single-particle wave functions and single-particle energies 
of an intrinsic state are calculated by diagonalizing the 
single-particle Hamiltonian matrix defined by analogy to 
the Hartree-Fock theory as done in Refs. nTolll7|. 

As 13 Be intrinsic states obtained with AMD+VAP are 
deformed and parity asymmetric, strictly speaking, spin 
and parity are not good quanta in each single-particle 
wave function. Nevertheless, single-particle wave func- 
tions can be associated with shell-model orbits from 
features of spatial distribution and ratios of positive- 
and negative-parity components. Figure 0] shows single- 
particle energies in the intrinsic wave functions of 13 Be 
states. For each single-particle level, fractions of posi- 
tive and negative-parity components are shown by a red- 
solid line and a grccn-dottcd one, respectively. Labels 



FIG. 4: (Color on-line) Single-particle energies in the intrinsic 
wave functions for 1J Be states, (a) 1/2", (b) 3/2", (c) 5/2", 
(d) 1/2+ , (e) 3/2+, (f) 5/2+, calculated with AMD+VAP us- 
ing the set (1) interaction. Fractions of positive- and negative- 
parity components are shown by red-solid and green-dotted 
lines, respectively. The labels, s, p, and sd indicate associa- 
tion with shell-model orbits. 



indicate rough correspondence to shell-model orbits. In 
the 13 Bc(l/2") state, the proton orbits and the lowest 
four neutron orbits form a la core. Among five va- 
lence neutrons around the 2cv core, three of them oc- 
cupy p-like orbits and the other two occupy sd-like or- 
bits. Thus, the 13 Bc(l/2~) roughly corresponds to a 
(p)~ 1 (sd) 2 configuration on the neutron p-shell and it is 
regarded asalfiw excited configuration. In a similar way, 
also the 13 Be(3/2~) and 13 Bc(5/2~) states have domi- 
nantly (p)~ 1 (sd) 2 configurations, and they are regarded 
as lhuj members as well as 13 Be(l/2~). On the other 
hand, the 13 Be(5/2 + ) contains mainly a (sd) 1 configura- 
tion and corresponds to the normal Ohuj configuration. 
In contrast to the 13 Be(5/2+), the 13 Bc(3/2 + ) contains 
dominantly a 2Tiu> excited configuration of a (p)~ 2 (sd) 3 
configuration. It is surprising that OHlu, ITiui, and 2Tiu> 
states almost degenerate in the low-energy region in 13 Be. 
For the 13 Be(l/2+) obtained with AMD+VAP, we can 
not associate valence neutron wave functions with shell 
model orbits because parities are strongly mixing and 
spatial behaviors show no specific feature analogous to 



shell-model orbits. Furthermore, this state may not cor- 
respond to the s-wave virtual state, which has been sug- 
gested near the 12 Be+n threshold. 
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FIG. 5: (Color on-line) Density distributions of the intrin- 
sic states for the Of and 0^ states of 12 Be calculated with 
AMD+VAP using the set (1) interaction. Same as Fig. [3] 
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FIG. 6: (Color on-line) Single-particle energies in the intrin- 
sic wave functions of 12 Be states. Same as Fig. 2] 

Let us compare deformations and cluster structures 
in 13 Be with those in 12 Be. Figures [5] and [6] show 
density distributions and single-particle energies of the 
intrinsic wave functions $( 12 Be; 2%u) for 12 Be(0 1 ") and 
$( 12 Be;0?iw) for 12 Be(0t) obtained with AMD+VAP. 
As discussed in Ref. [19|, 12 Be(0+) is dominantly the 
2Twj intruder state having a developed cluster structure 
with a large deformation, while 12 Be(0j) described dom- 
inantly by the normal Oftw configuration has a weaker 
cluster structure and a smaller deformation than those 
of 12 Be(0]|~ ). It is found that the large deformations 
in 13 Be(l/2-), 13 Be(3/2-), 13 Be(5/2"), and 13 Be(3/2+) 
are similar to that of 12 Be(0 1 "), while smaller deforma- 
tions in 13 Be(l/2 + ) and 13 Be(5/2 + ) are associated with 
12 Be(0^). Then, the 13 Be states can be interpreted as 
12 Be+n states by considering an aditional neutron on 
the 12 Be cores. That is, 13 Be(l/2"), 13 Bc(3/2-), and 
13 Be(5/2~) have the 12 Be(2huj) core with a neutron in 
a p-orbit. 13 Bc(3/2 + ) is regarded as the 12 Be(2?kj) core 
and a sd-orbit neutron, while 13 Be(5/2 + ) is interpreted 
as the 12 Be(0?kj) core and a sci-orbit neutron. Here cou- 
pling of the last neutron with the 12 Be(2?iw) core is not 
weak coupling but strong coupling where the neutron is 
moving around the largely deformed core. 



The 12 Bc+n features in the AMD+VAP results 
are supported also by the 12 Bc+n model calcu- 
lation. Indeed, in the 12 Be+n model calcu- 
lation, 13 Be(l/2~), 13 Be(3/2~), 13 Be(5/2~), and 
13 Be(3/2 + ) contain mainly 12 Be(2hui)+n components, 
while 13 Be(5/2 + ) and 13 Be(l/2 + ) are approximately de- 
scribed by 12 Be(0htu)+n wave functions. 



B. Molecular orbital picture 

In theoretical works on cluster structures of Be isotopes 
IHH, low-lying states of 9 Be, 10 Be, n Be, and 12 Bc have 
been successfully described in terms of molecular orbitals 
around a 2a core. As mentioned, the present work sug- 
gests the 2a core formation in 13 Be as well as in other 
Be isotopes. We here extend the molecular orbital de- 
scription to 13 Be and discuss structures of Be isotopes 
systematically. 

In the molecular orbital picture for a 2a system, molec- 
ular orbitals are formed by a linear combination of p 
orbits around two a clusters, and valence neutrons oc- 
cupy the molecular orbitals around the 2a core [8h13J|. 
A negative-parity orbital constructed by p-orbits per- 
pendicular to the a-a direction is called a "7r-orbital" , 
and a positive-parity orbital from p-orbits parallel to 
the a-a direction is called a "tr-orbital" (Figs. |7(a) and 
•ED 3 )). We call the other positive- party orbital given by 
p-orbits perpendicular to the a-a direction a "7r*-orbital" 
in analogy to electron orbitals in atomic molecular sys- 
tems (Fig.Efc)). 

To understand the breaking of magicity in a chain of 
Be isotopes, single-particle levels in the molecular orbital 
model (MO levels) is useful as discussed, for instance, in 
Ref. pjj]- In the molecular orbital models [9l-tl2l]. single- 
particle levels are evaluated as functions of the a-a dis- 
tance. They are smoothly connected from the one-center 
limit to the two-center limit, and those in the intermedi- 
ate region correspond to molecular orbitals. In addition 
to the spatial configurations (w, a, and it*), molecular 
orbitals are specified by the angular momentum O = j z 
projected on to the symmetric axis z. Q for the a- 
orbital is O = 1/2, and we use the notation <X\/2- For 
it- and 7r*-orbitals, O = 1/2 and 3/2 are possible. Due 
to the spin-orbit force, 7r-orbitals split into the Zs-favored 
7T3/2-orbital and the ^s-unfavored 7Ti/2-orbital, and it* 



orbitals do into the it 



3/2" 



and 7r*/ 2 -orbitals. 



'3/2 

Note that 



the present notations, 7r 3 / 2 , tt\ji, <J\[i, and 7r*, 2 corre- 
spond to the labels 7r3/2~(ff) a\/2~(g), <j1/2 + {u), and 
tt3/2+(m) in Fig. 15 of Ref. QJ], and the labels (3u, 1), 
(lu, 2), (lg, 2), and (lg, 2) in Ref. @, respectively. In the 
spherical shell model limit, the 7r 3 /2- and 7Ti/2-orbitals 
lead to the P3/2- and p^-orbits, respectively, while the 
CT1/2- and 7T3, 2 -orbitals go to the d 5 / 2 -orbit. The MO 
levels correspond well to the single-particle levels in the 
two-center shell model [4j|. Moreover, when a a-a dis- 
tance is not large, the MO levels are associated with the 



Nilsson levels of the deformed shell model [30]. 

One of the important features of molecular orbitals is 
that the cr-orbital has two nodes along the a-a direction 
and it gains the kinetic energy as the 2a cluster devel- 
ops. The breaking of N = 8 magicity in n Be and 12 Be 
is understood by the energy gain of the cr^-orbital in 
developed 2a systems. 

Let us describe structures of 13 Be with a 2a core 
and five valence neutrons in terms of molecular orbitals. 
In the developed cluster states of 13 Be obtained with 
AMD+VAP, neutron wave functions are associated with 
molecular orbitals. For instance, in the 1/2" state, three 
neutron wave functions correspond to the 7r-orbital and 
two neutron orbitals are associated with the cr-orbital. 
This state is described by a molecular configuration, 
7r 3/2 7r i/2 cr i/2' I n a s i mna r way, the 3/2 + state is un- 
derstood by a 7r 3/2 cr i/2 7r 3/2 configuration. In both cases, 
two neutrons occupy the <7i/2-orbital. Then, the reason 
for the low-lying intruder states in 13 Be can be under- 
stood again by the lowering cr^-orbital in the developed 
cluster structures. In other words, er^-orbital neutrons 
play an important role for the breaking of N — 8 magic- 
ity in 13 Be as well as 12 Be. As for the 5/2+ state in 13 Be, 
since it has a weaker cluster structure and a smaller de- 
formation, its neutron configurations should be associ- 
ated with spherical shell-model orbits rather than molec- 
ular orbitals. Nevertheless, taking into account the cor- 
respondence of the 773/2-, Tim-, and cr ]y 2 -orbitals to the 
P3/2-, Pi/2~i an d c?5/2-orbits, we here temporarily assign 

"-7T 2 -7T 2 

7r 3/2 7r l/2 

sion. 



7r i/2°i/2" t° the 5/2 + state in the following discus- 



Energy spectra of Be isotopes can be understood sys- 
tematically according to the ordering of MO levels. In 
the MO level ordering, the key feature is that such the 
orbitals as a- and 7r*-orbitals with nodal structures par- 
allel to the 2a direction gain kinetic energy in a devel- 
oped 2a system. As mentioned, the cr-orbital with two 
nodes gains the energy as the cluster develops. Also the 
7r*-orbitals with one node gain some kinetic energy in a 
developed cluster system. Differently from the a- and tt*- 
orbitals, the 7r-orbitals have no node and their energies 
increase relatively. Consequently, as the 2a cluster de- 
velops, the ordering of single-particle levels changes from 
spherical shell model orbits as shown in Fig. [7] (d). In 
this scenario, the breaking of the neutron magicity in 
Be isotopes occurs due to the intruder cr\ /2-orbital which 
comes down below the Tr^-orbital in the developed 2a 
systems. As a result, the N = 8 shell gap disappears 
and the level ordering based on a spherical shell model 
picture is no longer valid in neutron-rich Be. Instead, the 
level ordering of molecular orbitals works rather well to 
understand energy spectra of Be isotopes. From the fact 
that 11 Be(l/2 + ) and 11 Be(l/2~) almost degenerate, one 
may expect, in the first-order approximation, that the 
tti/2- an d <y 1 / 2 -orbitals almost degenerate in the MO lev- 
els. In the new ordering of the MO levels on a 2a core, the 
7r 3 / 2 -orbital should be the lowest, (J\i<x- and 7Ti/2-orbitals 



compose the second group (called a (Ji/2-^1/2 shell in 
the present paper), and the 7r*, 2 -orbital may exist above 
them (Fig. [3(d)). 

Let us review the molecular orbital configurations 
of Be isotopes. In Table IIVB1 configurations for va- 
lence neutrons around a 2a core for band-head states 
of Be isotopes are summarized. For 10 Be, lo Be(0^), 
10 Be(l~), and 10 Be(0 2 f ) are described by the valence 
neutron configurations of 7r 2 , 2 , 'K^i^xlii and tr 2 , 2 mean- 
ing two neutrons in the 7r-orbital, one neutron in 
the 7r 3 /2-orbital and the other neutron in the cr-- 
orbital, and two neutrons in the cr^-orbital, respec- 
tively. In a similar way, 11 Bc(l/2+), 11 Be(l/2-), and 
11 Be(3/2 2 ~) states are described by ^/ 2 < J i/2, 7r 3/2 7r i/2 
and ti'3/2 'i/2 configurations, while 12 Be(0f), 12 Be(0 2 f ), 
and 12 Be(lj~) states correspond to th^c 2 ^, 7r 3/2 7r i/2 
and 7t"3/2 7r i/2 cr i/2> respectively These assignments have 
been suggested by molecular orbital models [12|, \13^ , 
cluster models Jp, [H, HHH, and also the AMD 
model [H, [l7l4l9l |. For 13 Be, the molecular orbital 
configurations, 7 '' 2 /2 7r i/2 cr i/2 and 7r 3/2 cr i/2 7r 3/2> are as ~ 
signed to 13 Be(l/2~) and 13 Bc(3/2 + ), respectively, and 
7r 3/2 7r i/2 cr i/2 ^ s temporarily assigned to 13 Be(5/2 + ), as 
mentioned. For 14 Be, the ground state 14 Be(0+) is con- 
sidered to be a 7r2 /2 7r i/2' J i/2 configuration. T An ex- 
cited 0^ band was theoretically suggested by a VAP 
calculation [46[ and a 7r 3/2 (T i/2 7r 3/2 configuration was as- 
signed. 

Finally, we discuss energy spectra of Be isotopes in re- 
lation with molecular orbital configurations. In Fig. HI 
experimental excitation energies are shown for 10_12 Be, 
while theoretical values are shown for 13 Be and 14 Be[4q. 
In terms of spherical shell model levels, the inversion be- 
tween OTioj, lhu, and 2hco occurs in 11_13 Be, and the 
energy spectra seem to be out of the normal ordering. 
However, in terms of MO levels, the energy spectra of 
Be isotopes can be understood rather easily. In MO lev- 
els, the neutron Fermi level exists at the (J\/2-^i/2 shell 
in 10 ~ 14 Bc. Excited configurations are characterized by 
7r 3 / 2 -orbital holes or 7r*, 2 -orbital particles. Configura- 
tions without 7Tg / 2 particles nor 7r 3 / 2 holes are normal in 
the MO levels and they degenerate in a low-energy re- 
gion. For instance, the degeneracy of 11 Be(l/2 + ) and 
11 Be(l/2~) can be understood because all these states 
have "normal" configurations in MO levels. Also the 
coexistence of 12 Be(0+), 12 Be(0+), and 12 Bc(l-) in the 
low-energy region is to be expected because they have no 
excitation in MO configurations. In a similar way, the co- 
existence of 13 Be(l/2~) and 13 Be(5/2 + ) is not surprising 
as they have normal MO configurations. Excited config- 
uration states with 7^/2 holes or n* /2 particles have gen- 
erally higher excitation energies than low-lying normal 
MO configuration states. The numbers of 7r 3 / 2 holes and 
7Tg /2 particles are noted in Fig. [H] as well as Table. IIVBI 

10 Be(Oj), 10 Be(l"), 10 Be(0j), and u Bc(3/2-) have one 



or two 7r 3 /2 holes, and they exist in the excitation en- 
ergy region around E x = 4 ~ 6 MeV. The excited state 
with a 7Tg , 2 particle and that with two TTg , 2 particles are 



suggested in 13 Bc(3/2 + ) and 14 Bc(02 ), respectively. For 
more detailed discussion of energy spectra, one should 
take into account two-body correlations and a-a distance 
dependence of those MO levels. 
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FIG. 7: (Color on-line) (a)(b)(c) Schematic figures for molec- 
ular orbitals around a 2a core, (d) A schematic figure for 
evolution of single-particle level ordering from spherical shell- 
model levels to molecular orbital levels. 



TABLE I: Classification of band-head states in 10_14 Be. Har- 
monic oscillator shell-model configurations and molecular or- 
bital configurations are listed. The shell model configura- 
tions Ohcu, lhu, and 2Tiu are based on neutrons excited from 
the p-shell to sd-shell. For molecular orbital configurations, 
neutron configurations around a 2a core are described. The 
numbers of 713/2 holes and 71-3/2 particles are also shown. The 
molecular orbital configuration for 13 Be(5/2 + ) is a temporary 
assignment (See text). 
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FIG. 8: (Color on-line) Excitation energies of band-head 
states in 10_14 Be. Experimental data are shown for 10_12 Be, 
and theoretical values for Be[46[ are shown. The data for 
13 Be are the present results obtained by the 12 Be-n model 
calculation. Blue dotted, magenta dashed, and red solid lines 
indicate OTkJ, lhui, and 2Tiuj configuration states, respectively. 



Structure of 13 Be was investigated with VAP+AMD. 
In the AMD+VAP calculation using the set (1) inter- 
action that reproduces the parity inversion of n Be, an 
unnatural parity 1/2 - state was suggested to be lower 
than 5/2 + state indicating that the vanishing of N = 8 
magic number occurs in 13 Be. A low-lying 3/2 + state 
with a 2?io; configuration was also suggested. The present 
AMD+VAP calculation is a bound state approximation. 
To see effects of spatial extension of the last neutron wave 
function on the energy spectra, we also calculated 13 Be 
in the 12 Bc+n model. The degeneracy of the l/2~, 5/2 + , 
and 3/2 + states was suggested also by the 12 Bc+n model 
calculation. It is necessary to treat asymptotic behav- 
iors and out-going boundary more carefully to discuss 
detailed energy spectra of resonances. 

In analysis of intrinsic structures of 13 Be, large defor- 
mations with developed cluster structures are found in 
the intruder states such as the 1/2 - and 3/2 + states. 
These deformed states are regarded as states composed 
by one neutron on the deformed 12 Be core, which corre- 
sponds to the intrinsic state of 12 Bc(0f ). 

The intrinsic structures of 13 Be were also discussed 
in terms of molecular orbitals around a la core. The 



intruder states, 13 Be(l/2~) and 13 Be(3/2+) arc de- 
scribed by the molecular configurations, 7r 3/2 7r i/2 cr i/2 
and 7T3/2 (J i/2 7r 3/2' respectively. The breaking of N = 8 



magicity in 13 Be can be understood by molecular orbital 
levels as well as n Be and 12 Be. 
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